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ABSTRACT 

We have combined the results from our recent observations of Damped and sub- 
Damped Lyman-a systems with the MIKE and UVES spectrographs on the Magellan 
Clay and VLT Kueyen telescopes with ones from the literature to determine the Nh i- 
weighted mean metallicity of these systems based both on Fe, a depleted element in 
QSO absorbers and the local ISM, and Zn a relatively undepleted element. In each 
case, the Nh i weighted mean metallicity is higher and shows faster evolution in sub- 
DLAs than the classical DLA systems. Large grids of photoionisation models over 
the sub-DLA Nh i range with CLOUDY show that the ionisation corrections to the 
abundances are in general small, however the fraction of ionized H can be up to ~90 
per cent. The individual spectra have been shifted to the rest frame of the absorber 
and averaged together to determine the average properties of these systems at z < 1.5. 
We find that the average abundance pattern of the Sub-DLA systems is similar to the 
gas in the halo of the Milky Way, with an offset of ~0.3 dex in the overall metallicity. 
Both DLAs and Sub-DLAs show similar characteristics in their relative abundances 
patterns, although the DLAs have smaller ([Mn/Zn]} as well as higher ([Ti/Zn]) and 
([Cr/Zn]}. We calculate the contribution of sub-DLAs to the metal budget of the 
Universe, and find that the sub-DLA systems at z < 1.5 contain a comoving density 
of metals f2 OTe * ~(3.5-15.8)xl0 5 Mq Mpc -3 , at least twice the comoving density of 
metals in the DLA systems. The sub-DLAs do however track global chemical evolution 
models much more closely than do the DLAs, perhaps indicating that they are a less 
dust biased metallicity indicator of galaxies at high redshifts than the DLA systems. 
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1 INTRODUCTION 

Even with the recent introduction of large telescopes cou- 
pled with sensitive instrumentation, studying the properties 
of distant, high redshift (z) galaxies still remains a difficult 
task. When studied via emission, biases towards more lu- 
minous galaxies, i.e. the Malmquist bias, are often present. 
One of the most fundamental and basic predictions of galaxy 
evolution models is that the metallicity of the galaxy should 
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increase over time with the ongoing cycles of star formation 
and death that pollute the Interstellar Medium (ISM) of the 
galaxy with heavy elements in successively higher amounts 
with each subsequent generation of stars. 

Chemical evolution studies of the Milky Way are able 
to utilize individual stellar abundances over a wide range 
of metallicities to constrain models of the evolu t ion o f 
the Gala x y (see for instance | Pagel fc Tautvaisienel (|l99ct ); 
ISamlandl ij 19981 ); iNissen et all l|2000l . |2004| )). Direct stel- 
lar abundance measurements are difficult in even the most 
nearby galaxies, and currently impossible at higher redshifts. 
An opportunity to study the ISM of high redshift galaxies 
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occurs if a galaxy is fortuitously aligned along the line of 
sight to a distant Quasar (QSO), as the gas produces nar- 
row absorption lines in the continuum of the background 
QSO. 

These absorption line systems seen in the spectra of 
QSOs provide an alternative approach to studying distant 
galaxies, and offer a unique window into the ISM of high z 
galaxies as well a s the diffuse Interga la ctic Medium (IGM) 
(see for instance iLehner et alj (l2007h: iTripp et all 1 2008 ) ; 
iDanforth fe Shulll (|2008l l : iFrank. Mathur fc York! l|2008l ) for 
recent studies of the metal content of the IGM). As these 
galaxies are not selected upon any morphological criteria, 
this approach studies galaxies independent of a particular 
morphological type. The metal (Z>2) abundances in the 
ISM of these galaxies yield important clues into the pro- 
cesses of star formation and death and the overall chemical 
enrichment of the system being investigated. 

Based upon their neutral hydrogen column densities, in- 
tervening absorption line systems with strong Lyman-a lines 
showing damping wings can be classified into the Damped 
Lyman-a systems (DLAs with log Nh i ^2 0.3) and sub- 
DLA systems (sub-DLA 19 < log N H i < 20.3, |Peroux et all 
2001). The damping wings which give rise to the name of 
these systems and can be used to accurately measure Nh i 
in these systems begin the sub-DLA regime of log Nh i 
~19.0. With their high gas content, the DLA and sub- 
DLA systems are expected to be associated directly with 
galaxies at high redshift, and several syste ms have been 
confi r med with followup imaging at lower z dLe Brun et al.l 
Il997l : iBowen. Tripp fc Jenkinsll200ll : IChen et alJl2003h . 

After much investigation, the DLA systems have been 
found to have low metallicity (~0.15Zq) and show little 
chem i cal ev o lution from 0. 1 < z < 3.5 jKulkarni et al.l 
20051. 120071: iMeiring et all 120061 : IProchaska et all 120031 : 
Wolfe. Gawiser fc Prochaskall2005h 

tic chemical ev olution models (|F_ 

ICen et alJ Eool). If the DLAs, which do not appear to con- 
tain the metals predicted by chemical evolution models, 
could these missing metals be hiding in other locations, have 
they been systematically underestimated in DLAs, or both? 
This survey was conducted to obtain high resolution spectra 
in a homogenous sample of the lesser studied sub-DLAs to 
further investigate their role in cosmic chemical evolution. 

The global metal budget of the Universe can be es- 
timated based on the integrated comoving star formation 
rate density p» , and the theoretical m e tal production yield s 
of stars (p z ) (fSearle fc Sargentl fl972l ; IBouche et alj [2008). 
Adding the contributions from the Lyman-a forest, Ly- 
man break galaxies, su b-mm galaxies, BX/BM galaxies (i.e. 
lAdelberger et af]l2004h and the DLA systems to the global 
metal budget of the Universe, a dearth of ~ 40% of the ex- 
pected metal content is s een (IBouche. Lehnert fc Peroux I 
2006i; IBouche et al ] I2008I ). Several gro ups have reported 
high abundances in sub-DLA systems (iKhare et al J 120041 



. which contradicts galac- 
dPei. Fall fc Hauseri Il999l; 



Peroux et al. 2006a ; Prochaska et al. 2006; Me iring et alj 



20071 : ISrianand et al.ll2008l : iQuast. Reimers fc BaadefeoOSh . 
Strong Lyman-a forest lines with log Nh i <16 have also 
been recently observed to have super-solar metallicities, al- 
though ionisation c orrection factors are i mportant and not 
easily constrained l|Charlton et alj 120031 ; iDing et all 120031 ; 
M asiero et alj|2005l ). It has also been suggested that a large 



portion of these missing met als may be locked in to the hot 
intergalactic medium as well l|Ferrara et al.ll2005l ). 

Although redshifts 2 < 1.5 span 70% of the age of 
the Universe (assuming the "737" concordance cosmology 
of tt A = 0.7, Q. m = 0.3, H =70 km s" 1 Mpc" 1 which we 
utilize throughout this work), few observations have been 
made of z < 1.5 sub-DLAs due to the lack of spectrographs 
with enough sensitivity in short wavelengths and the the 
paucity of known sub-DLAs in this redshift range since Nh i 
must be determined with space based spectrographs as the 
Lyman-a line lies in the UV at these redshifts. This range is 
clearly important for understanding the nature of sub-DLA 
systems along with galactic chemical evolution. 

In this paper we discuss the full sample of sub-DLAs ob- 
tained in this survey with the Magellan-II telescope and the 
MIKE spectrograph and the VLT Kuyuen telescope with 
the UVES spectrograph. The structure of this paper is as 
follows. In §[2]we discuss details of our sample. In § 3, we dis- 
cuss the abundances and column density determinations. §2] 
gives information on photoionisation models with CLOUDY. 
In § 5 we discuss the kinematical properties of these systems. 
In §[5] we give the Nh i-weighted mean metallicity based on 
both Fe and Zn. We discuss the averaged sub-DLA spectrum 
in § [7] In § [8] we estimate the contribution of sub-DLA sys- 
tems to the comoving metal budget, and in §[5]we discuss 
the implications of these measurements. 



2 SAMPLE SELECTION 

The observations in this sample were made with the 6.5m 
Magellan-II Clay telescope and the Magell an Inamori Ky- 
ocera Echelle (MIKE) spectrograph (see iBernstein et al.l 
l2003h along with the 8.2m Very Large Telescope (VLT) 
Kueyen telescope wi th the Ultraviolet Visu al Echelle Spec- 
trograph (UVES, see lD'Odorico et al.|[2000l for details of the 
instrument) . These observations were undertaken from 2005 
until 2008. Our total sample includes spectra of 27 QSOs and 
a total of 34 systems, with 31 sub-DLAs or strong Lyman- 
limit sys tems and three DLAs. The s e resu lts h ave been pub- 
lished in lMeiring etafl (|2007l . 120081 . l2009h and lPeroux et all 

diooi). 

Our targets were selected from the 

iRao. Turnshek. fc Nestor] l|200et ) catalog of Mg II ab- 
sorbers identified in Hubble Space Telescope (HST) spectra 
having 0.65 < z abs < 1.5 (so that the Zn II AA 2026,2062 lines 
were redhsifted into the optical) with 19.0<Nh i <20.3. 
Due to the lower H I column densities and hence lower 
expected Zn II column densities for sub-solar metallicity 
systems in the Nh i range, we required the background 
QSO to be brighter than m g < 18.5 in order to achieve S/N 
ratios adequate to detect the expected weak Zn II lines. 

Ideally, a blind survey of QSOs is the best possi- 
ble option with minimal bias. This is however unfeasi- 
ble due to the necessity of both ground and space based 
spectra, and the small coincidence rate of Sub-DLAs or 
DLAs in the spectra of background QSOs. An Mg II se- 
lection criterion is not likely to significantly bias the sam- 
ple, as the rest frame equivalent width of the Mg II 
AA 2796, 2803 lines does not correl ate well with [Zn/H] 
or [Fe/H] for DLAs or sub-DLAs. dKulkarni et all [20071 ; 
iDessauges-Zavadskv. Ellison, fc Murphy||2009l ). Of the -45 
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Figure 1. A comparison of the rest frame equivalent widths of 
Mg II A 2796 and Fe II A 2 600 from these data, and from the 
lower resolution spectra from lRao. Turnshek. fc Nestorl (|2006h . 



systems from the iRao. Turnshek. fc Nestorl (l2006t ) sample, 
only eight systems remain at z<1.5 that fit the selection cri- 
teria above and not observed in the sample presented here 
due to scheduling and time constraints. These objects have 
similar mean Wq 796 (~1.4 A compared to ~1.6 A for sys- 
tems that we have observed), and thus are unlikely to bias 
the sample. 

Throughout this paper the QSO names are given in 
J2000 coordinates. We show a co mparison of the rest frame 
equiva lent widths as measured bv lRao. Turnshek. fc Nestorl 
|2006l ) and the values determined in these higher resolution 
spectra for the Mg II A 2796 and Fe II A 2600 lines in Figure 
[1] The values appear to be in excellent agreement, with no 
systematic differences in the two measurements. 



3 ABUNDANCES 

Total abundances for the observed systems are given in Ta- 
ble El These values have been given in our previous work s 
(see TMeiring et al.l 120071 120081 , I2OO9I ; iPeroux et all l2006bh , 
and are compiled here for completeness and convenience to 
the reader. We have used the total column densities (i.e. the 
sum of the column densities in the individual components 
of a system that were determined via profile fitting method) 
along with the total Nh i as given in Table [3] to determine 
the abundances of these systems. We have not placed any 
ionisation corrections while determining these abundances 
(the issue of ionisation corrections is discussed further in 
Section [4]), and have assumed the first ions to be the domi- 
nant ionisation species of the elements for which these abun- 
dances have been determined, nam ely Zn, Fe, Mn , Cr, and 
Si. Solar systems abundances from lLod dcrs (2003]) are also 
given in Table [3] 

In Table [3] we have given the metallicities of Zn and 
Fe, the relative abundances of a several species, as well as 
the column density ratios of adjacent ions of a few elements. 
All these are based on the column densities obtained from 
the profile fitting analysis. Relative abundances like [Zn/Fe], 
are indicators of dust depletion, and ion ratios like Al II/A1 



III, may provide information about the ionization in the 
systems. 



4 PHOTOIONISATION MODELS 

The gas in DLA systems is expected to remain largely neu- 
tral due to self-shielding of the UV photons capable of ion- 
ization. The ionisation correction factor, defined here as 



[X/H] t 



[X + /H°] 



(1) 



where the total column densities include contributions from 
all ionisation s tages, has been investigated for DLAs by 



all ionisation sta ges, has been mvesti g; 
several groups dHowk fc Sembach|[l999l ;T 
alT~ 



Vladilo et 
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IProchaska et all 2002). The ionisation corrections for most 



elements in DLA systems were found to be <0.2 dex in 
most cases. Due to the lesser amount of H I in the sub- 
DLA systems, and hence less shielding, it is expected that 
they should show a greater amount of ionisation. 

It has previously been shown that the ionisation correc- 
tions based on cloudy models a re in general also expected 
to be small for sub-DLA system s (iDessauges-Zavadskv et aU 
120031 ; iMeiring et al.ll2007l . |200S| ). Here, we expand upon our 
previous CLOUDY modeling with grids of photoionisation 
models over the entire sub-DLA Nh i range. Models were 
created with version C06.02.b of CLOUDY, last described by 
iFerland et al.l (l998). Grids of CLOUDY models were com- 
puted assuming that the spectrum of ionizing radiation 
striking the cloud followed the form of the extra-galactic UV 
background of Haardt fc Ma dau (|Haardt fc Ma dau 1996; 
iMadau. Haardt. fc R ccs 1999) at the appropriate redshift 
of the absorber, plus the model stellar atmosphere of Ku- 
rucz with a temperature of 30,000K to simulate a radia- 
tion field produced via an O/B-type star. Both radiation 
fields were combined in equal parts for the final incident 
spectrum. Plots of both types of spectra can be found in 
Hazy, the documentation for CLOUDY which can be found at 
\protect\vrulewidthOpi http://www.nublado.org . In addi- 
tion, our CLOUDY models have included a cosmic ray back- 
ground and cosmic microwave background. These simula- 
tions however do not include radiation from local shocks 
caused by supernovae, or compact sources such as white 
dwarfs or compact binary systems, all of which likely con- 
tribute the ionizing radiation field. 

For each of the grids of models, the ionisation parameter 
defined by 



(2) 



U=^L 

n H cn H 

(where $912 is the flux of radiation with hv > 13.6 eV) was 
increased from log U=— 6.0 to 0. 

Results for the ionisation corrections for Fe and Si are 
shown in Figures [2] and [3] as a function of Nh i and log U. 
As defined above, for a ratio X + / H° one must add e to get 
the true abundance X/H. Fe appears unaffected except in 
the extreme case of l ow Nh i and high ionisation param eter 
U. As was no t ed in Dessaug es-Zavadskv et al] l|2003l ) and 
IMeiring" et all l|2007l . |200S|) in all cases where a constraint 
on the ionization parameter could be placed in the modeled 
sub-DLA systems, log U< —2.5, indicating that the correc- 
tion factors are at most ~ —0.2 dex, in the sense that the 
measured abundances may be slightly higher than the true 
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Table 1: Elemental abundances and adjacent ion ratios for the absorbers in this sample. The solar abundances from Loddcrs (2003) arc given in the headings as well. In the last 
four columns, the ratios of the column densities of several pairs of ions are given in dex. Njj j is given in logarithmic form in units of cm -2 . 
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Figure 2. Results of the CLOUDY photoionisation models for the 
ionisation correction needed for Fe as a function of log U and 
Njj i- The blue contours on the horizontal plane represent the 
ionisation corrections in steps of 0.05 dex. 




Figure 3. Results of the CLOUDY photoionisation models for the 
ionisation correction needed for Si as a function of log U and 
Njj j. The blue contours on the horizontal plane represent the 
ionisation corrections in steps of 0.20 dex. 

abundances. With low ionisation parameters U, Fe + /Si + is 
an accurate indicator of Fe/Si, however at higher ionisation 
parameters (log U >-2), Fe + /Si + is higher than the true 
Fe/Si, mimicking the effects of differential and higher deple- 
tion of Fe than Si. 

We have also investigated the fractional ionisation of 
hydrogen in these sub-DLA systems to determine the over- 
all levels of ionisation which are needed to determine the 
contribution of sub-DLAs to the global metallicity. These 
simulations are more reliable than simulations of the ion- 
isation of the metal species due to the direct calculability 
of the atomic properties of hydrogen such as the recombina- 
tion rates, ionisation cross sections, and oscillator strengths. 
Again, we assumed a hot stellar spectrum with T = 30, 000 
K, and the diffuse extragalactic background from AGN and 
starburst galaxies at z ~ 1. Constant density was assumed 
with n = 1 atom cm -3 . Starting with log Nh i =18.8, the 
column density was increased by 0.2 dex in each subse- 
quent iteration. Three separate ionisation parameters were 




0.0 i i i i I I 

19 20 21 22 
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Figure 4. Results of the CLOUDY photoionisation models for the 
ionisation of hydrogen in these sub-DLA systems. The neutral 
fraction, H° /(H° + H + ) is plotted on the y axis, with log Njj i 
along the abscissa. Also plotted as red stars are the ionisation 
fractions derived from CLOUDY models of each individual system. 



used, with values of log U= —4.0, —3.5, and— 3.0. These val- 
ues were chosen because t h e sub - DLAs in this s a mple and 
iDessauees-Zavadskv et~aH l|2003h : IPeroux et all (|2007h all 
seem to have values of the ionisation parameter in or below 
this range. Results of the simulations are shown in Figure 
[4] Overplotted in the figure are the values for the ionisation 
parameters derived from individual CLOUDY models for the 
systems in this sample. 

In general, the sub-DLAs contain large fractions of ion- 
ized gas based upon these simulations, especially the systems 
with the lowest Nh i while the DLAs are largel y neutral in H. 
This was also reported for higher z systems by {Peroux et alj 
2007). However, for a given Nh i a large range of ionisation 
fractions are seen, perhaps indicating that there are several 
types of environments from which the sub-DLAs are being 
drawn. 

There are several unknowns however. The shape of the 
ionizing radiation field does affect the overall ionisation. The 
Haardt & Madau background spectrum is a theoretical con- 
struct, and has not been accurately tested and as such is 
only one of many possibilities. The hot stellar spectrum used 
in these simulations, although certainly measurable in the 
Milky Way, may also not accurately reflect the true radia- 
tion field incident on these systems. 



5 KINEMATICS OF SUB-DLAS 

Following the analysis of iLedoux et al.l {2006), we per- 
formed an analysis of th e apparent optical depth (see e.g., 
ISavage fe Sembachl [l996) of these sub-DLA systems. The 
apparent optical depth is defined as 

r(A)« =ln[I (X)/I obs (X)] (3) 

where Io(A) is the continuum level, and I f, s (A) is the ob- 
served intensity. Specifically, we integrated the AOD profile 
over the extent of the line, and give the velocity width (Augo) 
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as the value where the inner 90 percent of the absorption oc- 
curs. 

In Fig ure [5] the m etallic it y vs Avgp for the 



syste m s in ILedoux et all d2006h: iPeroux et al 



20081); iProchaska et al 



2006): Notcrdacme et al 



Quast. Reimers fe B aadc ( 20 08f) as well as the sys tems from 
this work. In their analysis, Ledoux et al.l (|2006h used ab- 
sorption features where 40 to 90 percent of the continuum 
level was absorbed by the line. We have tried to follow this 
strategy whenever possible for these systems. For most sys- 
tems, the Fe II A 2374 line was used in the velocity width de- 
termination, as it is typically unsaturated in these systems, 
and also strong enough to include absorption from weaker 
features. For systems where the Fe II A 2374 feature was too 
weak, the Fe I I A 23 44 feature was typically used. From the 
ILedoux et all l|2006h sample, we have only included the Zn, 
and not the Si and S detections in Figure [5] 

The non-parametric Wilcoxon Rank-Sum test was used 
to determine if the populations of sub-DLAs and DLAs were 
distinct in their mean velocity widths. The test statistic Z 
describing the significance of the difference in the means of 
the two populations was determined to be Z = 0.94, with a 
probability of the sample populations having the same mean 
of 17 per cent. Similarly, a KolmogorovSmirnov test shows 
that the populations of DLA and sub-DLA velocity widths 
having the same mean is ~18 per cent. The non-parametric 
approach used here is more appropriate than a two sample 
t-test, as we do not know the underlying distributions of the 
kinematics of th ese populations. A s imilar weak discrepancy 
was also seen bv lZwaan et al] (|2008l ) . wherein the lower col- 
umn density systems also had slightly larger Avqo values. 

As can be seen in Figure [5] the DLAs do show a fairly 
tight correlation between their metallicities and velocity 
widths. Interestingly, the sub-DLAs however do not show 
this trend. The Kendall's r test was used for both the DLAs 
and sub-DLAs to test for a correlation between the metal- 
licity and velocity width. For the DLAs, a clear correlation 
is present with r = 0.59 and a probability of no correlation 
p < 0.001. Conversely, the sub-DLAs show no evidence of 
correlation with r = 0.09 and a probability of no correlation 
of p ~ 60 percent. This however is possibly due to a lack of 
Zn detections at low [Zn/H]. 



6 N H i-WEIGHTED MEAN METALLICITY 
6.1 [Fe/H] Metallicity Evolution 

The Nh i-weighted mean metallicity given by: 
[(X/H)] = log((X/H)) - log(X/H) Q 
where 

Er.i 



with the errors estimated by the standard deviation, 



Y,M[X/H] t -[{X/H)]f 



/(« " 1) 



and the weights proportional to Nh i as given by 
_ N(H)j 



(4) 
(5) 

(6) 
(7) 
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Figure 5. The metallicity as measured by the undepleted ele- 
ments Zn and S vs the velocity width Avqq defined above. 



has often been used as a quantitative way of estimating the 
amount of metal enrichment in the Universe in a given epoch 
ijPeroux et al.ll2003bl : iKulkarni et aT]|2005l . l2007h . The N H i- 
Weighted Mean Metallicity is thus a measure of the global 
or average chemical enrichment of the Universe in a given 
epoch or redshift interval. 

As the lines of Fe are comparatively easier to detect 
and the cosmic abundance of Fe is substantially higher than 
that of the weaker Fe peak element Zn, Fe abundances for 
nearly every DLA and sub-DLA observed to date have been 
determined. There are also multiple lines of Fe II from ~ 
1200 to ~2600 A, allowing the possibility for higher redshift 
systems to be studied, where the Zn II AA 2026, 2062 lines 
would redshifted out of the visible range. On the other hand, 
as was previously mentioned, Fe is refractory and is typically 
depleted onto dust grains. [Fe/H] is therefore a lower limit 
to the true abundance. 

The Nh i-weighted mean metallicity based on Fe is 
shown in Figure [5] Sub-DLAs have been plotted as green 
squares, and DLAs as blue circles. These data were binned 
into five sub-samples of constant redshift size (Az ~ 1 
and solid points) and into bins with roughly equal num- 
bers of systems (open points). Independent of the binning 
type used, the sub-DLAs do appear to have a faster evolving 
metallicity, and are generally more metal rich in each red- 
shift interval. Al so shown in Figure [6] are the chemical evo- 
lution models of lPaeel fe Tautvaisiene] l l 19981) for the SMC 
(solid red line), and lPei. Fall, fe Hauserl 1 1999h (dashed ma- 
genta line) for the global metallicity evolution of the Uni- 
verse. Both types of models are s ophisticated open box ana - 
lytical models, with the models of lPei. Fall, fe Hauserl (|l999l ) 
being averaged over several comoving Mpc to contain mul- 
tiple galaxies. 



6.2 [Zn/H] Metallicity Evolution and Survival 
Analysis 

Although far easier to detect in QSO absorbers, Fe suffers 
as a metallicity indicator due to dust depletion. The non- 
refractory elements of O, S, and Zn are much more reliable 
indicators of the true metallicity of the system. The O and 
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Figure 6. The Nh i weighted mean metallicity for sub-DLAs 
(green squares) and DLAs (blue circles ) . O verlayed is the theoret - 
ical chemical evolution model from the lPei. Fall. fe Hauserl (|l999T ) 
model for the global metallicity (d ashed line), and the chemica l 
evolution model of the SMC from IPagel fe Tautvaisienel [|l998h 
(solid line). The lookback time in Gyr is given on the upper x 



S lines however lie far in the UV, at A ~1300 and 1250 A 
respectively and are inaccessible without space based obser- 
vations atz< 2.0. The Zn II AA 2026, 2062 lines however 
are accessible with ground based observations for these sys- 
tems. The nucleosynthetic origins of Zn are still somewhat a 
matter of debate. With Z=30, Zn is not quite a true Fe peak 
element. The exact production sit e of Zn is still in quest ion 
(see for instance the discussion in Mis henina et al.ll2002h . 

The Kaplan-Meier estimator has become the stan- 
dard non-parame tric statistical analysis for dealing with 
censored data l|Kaplan fe Meieil 1 1958h . The Kaplan- 
Meier estimator works with any underlying distribution 
of the data, which is advantageous in the case of this 
data set for which the underlying distribution of metal- 
licities at a given redshift or redshift interval is not 
known. We have used the ASURV package provided by 
the Astrostatistics Center at Penn State University 
(\protect\vrulewidthOpt http://astrostatistics.psu.edu ) 
for determining the summary statistics of the metallicities 
as measured by [Zn/H] for the DLAs and sub DLAs in this 
sample and the literature. The detail s of the package are 
described in lFeigelson fe Nelson! (|l985l ). for univariate data. 

The sub-DLA sample was separated into four bins with 
(roughly) equal numbers of systems in each bin, with the 
same done for the DLA systems from the literature along 
with the two DLAs from these observations. The DLAs were 
separated into five bins due to the larger sample size. A plot 
of the metallicity as determined via this survival analysis is 
shown in Figure [7] The vertical error bars denote the esti- 
mate of the variance in each redshift bin, and the horizontal 
bars denote the range of redshifts in each bin. Errors for the 
metallicity via survival analysis include only the statistical 
uncertainty and not any measurement uncertainties. As can 
be seen in both Figures [7] and [6] the metallicity based on 
Fe and Zn is faster evolving for the sub-DLA sample than 
the DLA sample. At z < 1.25, the sub-DLA metallicities are 
~ 0.7 dex higher than the DLA systems. 



Figure 7. The metallicity via survival analysis taking into ac- 
count upper limits for sub-DLAs (green squares) and DLAs (blue 
circles). Overlayed is the th e oretic al chemical evolution model 
from the |Pei, Fall, fe Hauserl jTJjlt) model for the global metal- 
licity (dash ed line) , and the chem ical evolution model of the SMC 
from Pagel & Tautvaisienc (1998) (solid line). The lookback time 
in Gyr is given on the upper x axis 



7 AVERAGE SUB-DLA SPECTRUM 



lYork et al.l (|2006h used a sample of ~ 800 lower resolution 
(~2000) SDSS spectra to determine abundances and deple- 
tion levels in QSO absorbers by averaging the spectra after 
shifting them to the rest frame of the absorbers. Weaker 
absorptions, not seen or seen at lower significance in the 
individual SDSS spectra, were clearly detected in this av- 
erage spectrum owing to its higher S/N ratio. Following in 
the spirit of that work, we have also averaged the spectra of 
these sub-DLAs. 

To determine average abundances of Cr, Fe, Mn,Si, and 
Ti of the absorbers in this sample, the spectra were shifted 
to rest wavelengths and averaged. Specifically, the Si II A 
1808, Cr II A 2056, Fe II AA 2260, 2374 and Zn II A 2026 
lines were used. Any obvious non-associated features in the 
region around the lines of interest in the individual spectra 
were replaced by Gaussian noise with S/N similar to the un- 
contaminated regions. For Q0251+1950, a systematic shift 
of ~60 k m s~ was applied to the spectra as the redshift 
reported in lRao. Turnshek. fe Nestorl (|2006l ) was slightly off. 
All other spectra were left unaltered. We show the averaged 
spectra of these transitions in Figure [8] 

Column densities for these elements were determined 
via the AOD method. The profile fitting method, although 
more robust and generally more accurate can not be used in 
this the individual components are blended together 

in the averaged spectra. Rest frame equivalent widths were 
also determined, with the errors estimated by both photon 
noise and uncertainties in the continuum placement. Abso- 
lute abundances were determined by averaging Nh i for the 
systems that were included in the averaged metal line spec- 
trum, while excluding the systems which were not used in 
the average. Table[5]gives Wo, the measured column density 
and corresponding abundance. Interestingly, the metallicity 
based on Zn from this averaged spectrum matches nearly 
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Figure 8. The average spectra of the systems in this sample. 
The key lines of Si II, Fe II, Zn II, Cr II, Mn II are shown here. 



Table 2. Equivalent widths, abundances and column densities of 
the averaged spectra. 



Transition 


# Systems 


Wo 


log N 


[X/H] 






mA 


cm 


Dex 


Si II 1808 


15 


44±9 


14.87±0.08 


-0.52 ±0.08 


Ti II 3383 


15 


12±4 


11.52±0.17 


-1.14 ±0.17 


Cr II 2056 


29 


15±4 


12.59±0.12 


-0.91 ± 0.12 


Mn II 2576 


25 


51±8 


12.39±0.06 


-0.98 ± 0.06 


Fc II 2260 


29 


37±5 


14.53±0.06 


-0.79 ± 0.06 


Fc II 2374 


28 


327±7 


14.38±0.01 


-0.93 ± 0.03 


Zn II 2026 


31 


28±5 


12.10±0.08 a 


-0.38 ± 0.08 a 



"Based on the average column density determined with and with- 
out the feature at ~ ±60 km s — 1 , which could be due to a con- 
tribution from Mg I 2026. 



perfectly with the Nh i-weighted mean metallicity from Sj6l 
validating the technique. 

The abundance patterns seen in the ISM of the Milky 
Way, SMC, and DLAs are shown i n Figure [9] Abund ances 
from the ISM of the SMC are fro mlWeltv et all dl997h wit h 
abundances from the DLAs from (|Prochaska fc W olfe 2002). 
Several points come to light based on the depletion patterns 
and abundances of these samples. The ISM Si abundance is 
nearly identical in all three cases. Secondly, the Zn abun- 
dance of the averaged sub-DLA spectra lies between the 
SMC and warm Galactic ISM value. Lastly, the Fe and Cr 
abundances in the averaged sub-DLA spectra are somewhat 
higher (~ ±0.5 and ~ ±0.3 dex for Fe and Cr respectively) 
than what is observed in the warm neutral ISM of the Milky 
Way and SMC. 

As the abundance patterns seen in the ISM are closely 
related to the past star formation rate and history of the 
galaxy, does the average abundance pattern in sub-DLAs 
mimic the patterns seen in the Milky Way or ISM, or nei- 
ther? The SMC and Milky Way have had very different star 
formation histories, and thus have different elemental abun- 
dances. The star for mation history of the SM C has recently 
been investigated bv lHarris fc Zaritskvl (|2004l ) based on deep 



photometric surveys of the entire SMC. They find that most 
of the stars in the SMC formed >8.5 Gyr ago (corresponding 
to redshift z > 1.2) with a period of quiescence lasting ~3 
Gyr, followed by more recent bursts in star formation. Tidal 
interactions with the Milky Way li kely play a roll in causin g 
star formation bursts in the SMC (|Harris fc Zaritskvll2004 ). 
The Milky Way has also gone through bursts of star forma- 
tion, creating the halo and thick disk. The Milky Way how- 
ever does not appear to have undergone the same long dura- 
tion periods of quiescence that the SMC und erwent, but in- 
stead has seemed to con tinuously form stars l|Gilmoreil200ll ; 
iRocah-Pinto et aI1l200Qt ). Sub-DLAs may trace a population 
of galaxies dissimilar to massive spirals (i.e. the Milky Way) 
or lower mass irregulars (i.e. the SMC) and thus may have 
different abundance patterns. 

Based upon the average abundance pattern in these 
sub-DLA systems, is appears as though these systems are 
more closely related to the Milky Way halo type interstel- 
lar medium than the warm neutral ISM of the Milky Way. 
The Zn metallicity in sub-DLAs, [Zn/H], as measured in the 
averaged spectra lies between the values seen in the warm 
neutral ISM and SMC. [Zn/H] is ~0.2 dex less than that of 
the warm Milky Way ISM and 0.3 dex higher than the SMC. 
This however reflects [Zn/H] at z ~ 1.1, the average redshift 
of these systems. This metallicity is well within the range of 
values at this redshift from the chemical evolution models 
mentioned above. The abundance and depletion pattern of 
the sub-DLAs does very closely resemble that of the halo of 
the Milky Way, as can be seen in Figure [lOl 



8 THE CONTRIBUTION OF SUB-DLAS TO 
THE METAL BUDGET OF THE UNIVERSE 

It has now been well established that the DLA systems 
observed to date are metal poor at all redshifts, and 
show weak evolution w i th redshift ([Kulkarni et al.l 120051 ; 
iProchaska fc Wolfdl2002l ; lMeirmg et al.ll2006T ). Including the 
contributions from star forming Lyman-break galaxies at 
2 ~ 2, DLAs, and the Lyman-Q forest the observations ac- 
count for only ~ 60% of the amount of metal s that are ex- 
pecte d from the cosmic star formation history l|Bouche et al.l 
2008). Here, we discuss the contribution of the less studied 
sub-DLA systems to the metal budget of the Universe. For 
a previous estimation of the comoving metal density in sub- 
DLAs based on the preliminary res ults of this work and t he 
literature sample of sub-DLAs, see [Kulkarni et all ([2007) . 

We can estimate the comoving metal density in DLAs 
and sub-DLAs based on the column density distribution 
function, metallicities, and in the case of sub-DLAs the ion- 
isation fraction. As was discussed previously, and can be 
seen in Figure [4] the ionisation fraction in the DLA systems 
with log Nh i >20.3 is negligible, and they can be considered 
completely neutral to good approximatio n in the H I phase. 
See however I Wolfe and Prochaskal (|2000h for a discussion on 
the more highly ionized gas in DLAs. 

The DLAs are well sampled, especially at z > 1.5, and 
have a mean metal licity of [Zn/H] ~-1.0 from a robust sample 
of ~ 1 00 systems (IProchaska fc Wolfel |2002| ; iKulkarni et all 
120051 ; iMeiring et alj|2006ri. Using the solar metallicity by 
mass as determined by Lodders! (2003) of Zq=0.0177, the 
DLAs thus have a metallicity by mass of Zdla=0. 00177. 



A Sample of Sub-DLAs at z < 1.5. 



9 



3olar_(100% Gas_ Phase) 




-2.0 - 



Zn Si Mn Cr Fe Ti 

Element 

Figure 9. The depletion pattern or abundances derived from the 
averaged sub-DLA spectra of the systems in this sample. Over- 
plotted are the abundances of the Milky Way warm neutra ISM, 
Milky Way halo, ISM abundances in the SMC, and DLA abun- 
dances. 
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Figure 10. The abundances relative to Zn derived from the 
averaged sub-DLA spectra of the systems in this sample. Over- 
plotted are the abundances of the Milky Way warm ISM, Milky 
Way halo, ISM abundances in the SMC, and DLA abundances. 

The fraction of mass in metal ions in DLA systems compared 
to the critical density of the Universe is thus 

c\DLA Pmet ^.DLA ~ / Q \ 

"met = = S1 H7 1 ^DLA (,») 

pcrit 

iPeroux et all (|2005h measured fifj/' 4 at z ~ 2.5, where the 
bulk of DLA metallicities are determined, to be fi^/* 4 ~ 
0.85 x 10" 3 . This then gives the comoving mass density of 
metals in DLAs to be O^et 4 ~2x 10 5 M n Mpc" 3 (see also 
iKulkarni et~al]|2007l ; IProchaska et al.ll2006h . 

For sub-DLA systems, the ionisation fraction of the gas 
must be taken into account to determine the comoving mass 
density of metals contained in these systems. If we let y(Nm) 
be the neutral fraction of gas in sub-DLAs (which is a func- 
tion of the neutral hydrogen column density as can be seen 
in Figure 2]), f(Nm) the column density distribution in the 
sub-DLA regime, and the average metallicity Z su tDLA we 
have then, 



n subDLA = M#omH f ± €{Z) f{Nm) z subNm dN m (9) 

Cpcrit J V 

where e(Z) is the ionisation correction term, which as can be 
seen from Figures [2] and [3] is a function of Nh i and the inte- 
gral runs over the sub-DLA column density range of 19.0 < 
log Nh i < 20.3. We have used the ionization correction de- 
termined for Fe for the purpose of these calculations due to 
the better quality of atomic data for Fe. 

Unfortunately, at lower redshifts the comoving gas den- 
sity in sub-DLA systems is not well known, owing to the lack 
of sp ace based UV survey s for H I at low z. Recently h ow- 
ever. lO'Meara et al.l (|2007h and also lPeroux et al.l l|2005h de- 
termined the column density distribution function (CDDF) 
and comoving mass density of gas in sub-DLAs from obser- 
vations with the Keck ESI, Magellan II MIKE spectrograph, 
and the VLT Kueyen UVES spectrograph. Combining the 
observations from both teams, we find that the CDDF is 
well fit by a power law (f{N m ) = with [3 = 10 4 84 

and a — —1.30 in the sub-DLA regime of 19.0 < log Nh i < 
20.3. 

Survival analysis provides an estimate of the Nh i- 
weighted mean metallicities given in Table [3] Using these 
values for the column density distribution function and mean 
metallicities, Equation [9] can be integrated. As can be seen 
in Figures [6] and [7] the sub-DLA and DLA populations are 
dissimilar in their respective metallicities, sub-DLAs hav- 
ing a significantly higher Nh i-weighted mean metallicity of 
{[Zn/R] subDLA ) ~ -0.08 and {[Zn/B] DLA ) ~ -0.8 at z < 1. 
These calculations were carried out using three ionisation 
parameters and thus three separate levels of ionisation for 
these systems. We have used the results from the photoioni- 
sation models for Fe for the correction term ^(A^HOFeii, when 
evaluating Equation [9] 

Results for the calculations of the comoving mass den- 
sity of metals in the sub-DLAs are given in Table [3] The 
sub-DLAs thus appear to contain several times the comov- 
ing mass density of metals that the DLA systems contain at 
z < 1, and at least equal amounts at all redshifts. We stress 
however that the comoving mass density of H, Q^ A , is not 
well known at low z and warrants further investigation. 

The more highly ionized phase of sub-DLAs has also 
been investigated as a possible location of the missing met- 
als. At low density and temperatures T ~ 10 6 K, the cooling 
time of the gas is roughly the age of the universe , so metal 
atoms can be essentially stuck in this hot phase. iFox et all 
( 2007) estimates that the metal content of the hot plasma in 
sub-DLAs at z ~ 2.5 is roughly equal to the what we have 
determined here for the lower temperature gas in which Zn 
II arises. 



9 DISCUSSION AND CONCLUSIONS 

We have compiled all of the known Zn measurements in 
DLAs and sub-DLAs, and have examined log Nh i vs. 
[Zn/H]. As has been reported in the past, the observed 
trend towards increasing metallicity with decre asing Nh i is 
again seen with the inclusion of these new da ta (|Boisse et al.l 
ll998l ; lKhare et al.ll2007l ; llvleiring et al.ll2006l ) . We have used 
survival analysis to calculate a trend line to these data 
including the upper limits, which is determined to be of 
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Rcdhsift 


Metallicity 




Pz 


Range 


Dex 




Mq Mpc -3 


0.53 - 0.87 


+0.06 ±0.22 


(7.8-31.3) xlO" 6 


(10.6-42.5) xlO 5 


0.89 - 1.14 


-0.29 ±0.17 


(3.5-13.9)xl0~ 6 


(4.7-18.9) xlO 5 


1.15 - 1.55 


-0.68 ±0.20 


(1.4-5.7)xl0~ 6 


(1.9-7.8) X 10 s 


1.77 - 3.39 


-0.59 ±0.22 


(1.8-7.0)xl0- 6 


(2.5-9.5) xlO 5 



and 



Table 3. The comoving mass density of metals, 
pz in sub-DLAs. 

a Calculated based on the z ~ 2 column density distribution 
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Figu re 11. Log Nh i vs. [Zn/H] for all the systems given in 
iKulkarn i et al.l 1120071), systems from this sample, and system s 
fromlQu astT^eimers fe Baadd d2008h;l Noterdaeme et all (l2008al) ; 
iProchaska et al.1 [|2006h ; iDessauges-Zavadskv et al.l d2003h. The 
solid line is the "Obscuration Threshold" from iBoisse et al.l 
(1998), while the dashed line is the linear fit to these data ig- 
noring the [Zn/H] upper limits, and the dashed-dotted line is the 
linear fit while including the upper limits. 



the form [Zn/H]= -0.58N H i + 10.9110. Without inclu- 
sion of the upper limits in the trend line, we determined 
[Zn/H]= -0.78N H i ±15.27. Both Spearman and Kendall r 
tests independently determined that the probability of no 
correlation is p < 0.001, with t he inclusion of th e upper lim- 
its. A similar trend was seen in lYork et al.1 (|2006l ) using com- 
posite spectra of a large number of SDSS QSO absorbers, 
and in iKhare et al l (2007) using the literature data as well. 

Interestingly, the amount of dust as estimated by 
[Fe/Zn] seems independent of Nh i in QSO absorbers over 
~3 decades in Nh i as is shown in Figure 1121 Spearman 
and Kendall correlation tests each reveal a probability of 
no correlation of p ~ 84 percent. The dust depletion on 
the other hand is known to correlate with the metallicity of 
the system as shown in Figure [TS1 where we show [Fe/Zn] 
vs [Z n/H] for QSO absorbers with both Fe and Zn detec- 
tions (|Kulkarni et al. 2007 ; Dessauges-Zavadskv et al. 2003; 



iQuast. Reimers fc Baadd (20081 ; iNoterdaeme et al.l l2008aF ) 
Spearman and Kendall correlation tests show that the prob- 
ability of no correlation between the metallicity and deple- 
tion as estimated by [Fe/Zn] to be p < 0.001. 

The observed trend of increasing metallicity with de- 
creasing Nh i seen in Fig. [11] contradicts measurements in 
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Figure 12. [Fe/Zn] vs. Log Nh i for QSO absorbers with 
both Fe an d Zn detection s from t his work, as well as sys- 
tems from [ K ulkarni et al.l l l2007l); Dessauges-Zavadsk v et all 
i2003lh-|Prochaska et al.1 [ 20061) ; IQuast. Reimers fe Baadd i200St) ; 
Noterdaeme et al.l ll2008al) ! 
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the local ISM of the Milky Way and the Magellanic Clouds, 
where different lines of sight show roughly solar metallicity 
gas although the sightlines may contain a range of Nh i val- 
ues. Abundance gradients could not produce the effect as 
is seen; if the majority of DLAs are produced in the out- 
skirts of the disks of massive spirals where the abundances 
would be lower due to the gradient, a simultaneous gradient 
of decreasing Nh i with decreasing galactic radius which is 
not observed would be necessary to explain the lower Nh i 
systems with higher abundances. 

In all likelihood, the sightlines through DLAs and sub- 
DLAs pass through gas rich galaxies of a wide range of 
morphologies and masses. The fact that the ionization frac- 
tions in sub-DLAs from the models above as well as in 
IPeroux et aD (|2007n show a range of values for a given Nh i, 
as well as the fact that not all sub-DLAs are metal rich sug- 



gest that this is indeed the case. There are however several 
lines of evidence to suggest that DLAs are mainly sampling 
low mass gas rich dwarf galaxies such as the SMC: 

• The abundances in DLA systems are low ([Zn/H]~ 
—0.8) and slowly evolving, which is consistent with the 
observations and chemical evolution models of local dwarf 
galaxies. The Nh i-weighted mean metallicity of DLAs does 
not track the global chemical evolution models and is es- 
sentially flat for 1 < z < 3, which is predicted by 

iPagel fc Tautvaisiend (|l998l ) for the SMC. 

• The molecular fraction of gas in DLAs is small 
iLedoux. Petitiean. fc Srianandl ll2003h. consistent wi th CO 
surveys of local dwarf galaxies ( Lerov et al.l 1200 5) . The 
molecular fraction in dwarf ga laxies is substanti ally smaller 
than in more massive systems (jLerov et al.ll20 05). consistent 
with FUV stellar obse rvations in the Magellanic Clouds (e.g. 
iTumlinson et al .120021 ). Although H 2 was detected in most of 
the SMC sight lines from lTumlinson et all i|2002l ). the molec- 
ular content of the SMC was likely much lower in the past as 
there was a long quiescent period in the star formation his- 
tory of the SMC. The molecular content of sub-DLAs has not 
been investigated as thorou ghly, however molecule rich sub- 
DLAs have been observe d (|Quast. Reimers fc Baadel 120081 ; 
iNoterdaeme et al.|[2008al '). 

• Followup imaging studies of DLAs show that the major- 
ity of the galaxies assoc iated with the DLA absorption, when 
detected have L < L * (|Le Brun et al.lll997l ; rRao et alj|2003l ; 
iKulkarni et al.ll2006l'). The lo w mass local dwarf galaxies ob- 
served in iGeha et alj (|2006l ) were determined to have total 
baryonic masses actually dominated by the gas, and not the 
stellar content, with gas fractions as high as f gas — 95% 
and {fgas} ~ 0.6. This could explain the low luminosities 
observed in DLA galaxies, and the low detection rate in fol- 
lowup imaging. 

This still raises the question as to why we do not see 
metal rich DLAs and do see more metal rich systems at 
lower Nh i, and the observed relationship between [Zn/H] 
and Nh i. One proposed explanation is that there is a selec- 
tion effect causing the higher metallicity DLAs to be system- 
atically undersampled due to dust obscuring the background 
QSO to magnitudes fainter than wha t is observable with the 
current classes of instrumentation ( Lauroesch et alj 1 19961 ; 
iBoisse et alll99Sl ; IVladilo fc PerouxH2005l) . The amounts of 
dust in QSO absorbers as estimated by ratios of volatile and 
refractory elements such as [Cr/Zn] or [Fe/Zn] is strongly 
correlated with the systems metallicity (Me iring et al.|[2006l ) 
the higher metallicity systems thus contain more dust and 
could be obscuring the background QSO. 

The amount of extinction need not be very large to 
exclude a large number of background QSOs from being ob- 
served. In fact, as can be seen in Figure [T4l only ~ 4 percent 
of SDSS DR7 QSOs have m s , < 18.0 at 1.25 < z em < 3.5, 
so even with a modest 0.5 magnitudes of extinction (less 
than 50 percent o f the estimated extinction amounts in 
IVladilo fc Perouxl {2005) for a solar metallicity DLA) we 
would be left with only 4 percent of QSOs bright enough to 
observe. As we are currently magnitude limited at m g < 18.5 
with current instrumentation, a large number of DLAs could 
be lying unobserved in front of fainter QSOs. 

If this is indeed true, then one would expect that the 
observed DLAs would have properties more similar to dwarf 
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Figure 14. The distribution of SDSS g' and r' magnitudes of 
QSOs with 1.25 < z < 3.5. The vertical lines at m =18.0 and 
m =18.5 represent the approximate magnitude limitations of the 
current instrumentation with and without an extinction effect 
respectively. 



or SMC type galaxies. The DLAs do seem to exhibit a more 
gradual chemical evolution similar to the SMC, as seen in 
Figures [6] and [7] and the molecular content along with the 
low luminosities of DLA host galaxies observed in imaging 
support the claim, as highlighted in the list above. With 
lower Nh i and thus likely lower levels of extinction, the 
sub-DLAs are less affected by the d ust extinction bias tha n 
the more H I rich traditional DLAs (jLauroesch et al.1ll996l ). 
This could possibly explain why the mean abundances in 
these sub-DLA systems are significantly higher than their 
DLA count erparts. An a l ternat ive explanation has been put 
forward bv lKhare et alj l|2007h . who sugges t that based on 
the ob served mass metallicity relationship of lTremonti et ail 
(2004) these higher metallicity sub-DLAs may arise in more 
massive spiral galaxies while the lower metallicity DLAs 
would be more similar to gas rich dwarfs. 

More work still is needed in several areas though. More 
space based observations are necessary to constrain the low 
redshift column density distribution of sub-DLA systems, 
and hence Qf^j bDLA which is currently poorly understood. If 
QsubDLA var j ec i f rom z ^ 2.5 to z ~ 1.0, the contribution 
of these sub-DLAs to the metal budget will alter accordingly. 

The Lyman limit systems with 17.2<log Nh i 
<19.0 are yet another possible reservoir of metals 
that have been hitherto understudi ed, although p r ogress 
is being made (see for instance [Charlton et alj 120031 ; 
M isawa. Charlton fc Narayanan! I2008T ). Several difficulties 
arise in these systems as Nh i cannot be constrained due 
to the lack of damping wings in the Lyman-cv profile in this 
column density range, and with such low neutral hydrogen 
column densities the effects of ionization are likely to be 
much more important. Nonetheless, this is certainly an area 
worth investigation. 

The ionisation levels in sub-DLAs are still not particu- 
larly well known, which is why we have presented a broad 
range of cases to illustrate the range of possible values in 
determining fi me t and pz- With further UV spectra of sub- 
DLAs, the higher ionization lines of S III, Si III, Fe III could 
be accessed, providing further constraints on the levels of 
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ionization in the gas. The Si Ill/Si II ratio seems to be a 
promising ionization indicator, as the ratio increases faster 
than other adjacent ion ratios such as Fe III/Fe II, enabling 
tighter constraints to be place on the ionization parameter. 
Simultaneously, the H2 transitions could also be observed to 
determine the molecular content in these systems. 

This series of papers has roughly doubled the total num- 
ber of sub-DLA observations, and nearly quadrupled the 
sample at z < 1.5, and done so with quality high resolution 
spectra. Complementary to studying these galaxies in ab- 
sorption, the typically more difficult task of imaging these 
galaxies still needs to be undertaken. Deep imaging in mul- 
tiple bands (g', r', i', K) of sub-DLAs is necessary to deter- 
mine the stellar content of these systems. Almost certainly, 
the sightlines that produce sub-DLAs are sampling a variety 
of morphological types of gas rich galaxies. The systems with 
high metallicity and higher velocity widths may however be 
more representative of massive spiral galaxies. 
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